Diffusion-controlled annihilation A + B — > with initially separated reactants: The 
death of an A particle island in the B particle sea 
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We consider the diffusion-controlled annihilation dynamics A+B — » with equal species diffusivities 
in the system where an island of particles A is surrounded by the uniform sea of particles B. We 
show that once the initial number of particles in the island is large enough, then at any system's 
dimensionality d the death of the majority of particles occurs in the universal scaling regime within 
which m 4/5 of the particles die at the island expansion stage and the remaining « 1/5 at the stage 
of its subsequent contraction. In the quasistatic approximation the scaling of the reaction zone has 
been obtained for the cases of mean-field (d > d c ) and fluctuation (d < d c ) dynamics of the front. 
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The reaction front formed in the A + B — > reaction- 
diffusion systems with initially separated reactants is of 
great interest since it represents a pattern for a wide spec- 
trum of processes in physics, chemistry, and biology [1]. 
Since the seminal article of Galfi and Racz [2] much work 
has been devoted to studying this problem by different 
approaches [3-20]. A standard way to treat the problem 
analytically is to solve the system of equations 

da/dt = D A V 2 a- R, db/dt = D B V 2 b - R (1) 



with the initial state given by 

a(x, 0) = clq6(— x), b(x, 0) = bo6(x), 



(2) 



where a(x,t) and b(x,t) are the mean local concentra- 
tions of A' s and B's, R{x,t) is the macroscopic reaction 
rate and 6(x) denotes the Heaviside step function, so that 
the A's are initially uniformly distributed on the left side 
(x < 0), and the B's on the right side (x > 0) of the ini- 
tial boundary. 

Dimensional [9,14] and renormalization group [11,12] 
analyses show that at d > d c — 2 one can adopt the 
mean-field approximation R(x,t) — ka(x,t)b(x,t)(k be- 
ing the reaction rate constant) with logarithmic correc- 
tions in the 2d case, whereas in Id systems fluctuations 
play the dominant role and the explicit form of R remains 
unknown. There are, however, several techniques which 
enable one to derive a lot of information from (1) even for 
d < d c . They are focused on the long-time limit kt — > oo 
and include, as two basic concepts, the scaling ansatz [2] 
and the quasistatic approximation [8,9,11]. According to 
the scaling ansatz (SA) the long-time behavior of the sys- 
tem inside the reaction zone may be represented in the 
form 



R = R/Q 



(3) 



where Xf cx t 1 / 2 denotes the position of the reaction zone 
center, Rj cx t^ 13 is the height, and w cx t a is the width of 
the reaction zone. At d > d c the scaling exponents take 
the values a = 1/6,(3 = 2/3 [2] whereas at d = 1 they 



appear equal to a = 1/4, (3 = 3/4 [9-14], so that at any 
d the width of the reaction zone on the diffusion length 
scale Ajj cx t 1 / 2 asymptotically unlimitedly contracts: 

w/ K-d — > as t — * oo. 

The quasistatic approximation (QSA) consists in the as- 
sumption that for sufficiently long times the kinetics of 
the front is governed by two characteristic time scales. 
One time scale, tj — —(dlnJ/dt)^ 1 , controls the rate of 
change in the diffusive current J = J a = Jb of parti- 
cles arriving at the reaction zone. The second time scale, 
tf oc w 2 1 D is the equilibration time of the reaction front. 
Assuming tj/tj <C 1 from the QSA in the mean-field case 
with D a ,b = D it follows [8,9] 



R 



I 



J/w, w - (D 2 /Jk) 1 / 3 , 



(4) 



whereas in the lei case w acquires the fc-independent form 
w ~ (D/J) 1 / 2 [9,11,17]. The most important feature of 
the QSA is that w and R / depend on t only through the 
time dependent boundary current J(t), which can be cal- 
culated analytically without knowing the concrete form 
of Q, i.e., in fact, representing the reaction zone on the 
scale Arj in the form R(x,t) = JS(x — Xr). On the ba- 
sis of the QSA the general description of spatiotemporal 
behavior of the system A + B — > has been obtained for 
arbitrary nonzero diffusion coefficients and initial species 
concentrations [19]. These results are in full agreement 
with extended numerical calculations and experiments 
and were generalized recently to the cases of reversible 
reaction A + B <-> C [21-24] and to several more complex 
reactions [25-29]. 

Until now, however, the main attention has been fo- 
cused on the systems with A and B domains unlimited 
extension, i.e., with unlimited number of A's and B's 
particles. The aim of this work is to develop a new line 
in the study of the A + B — > dynamics under the as- 
sumption that the particle number of one of the species, 
say A, per unit of the initial boundary is finite. More pre- 
cisely, we will consider the problem on the dynamics of 
death of an A particle island surrounded by the uniform 



1 



sea of particles B and will reveal the denning features of 
this process. 

Let particles A with concentration cio be uniformly dis- 
tributed in the island x € (-L, L) surrounded by the un- 
limited sea of particles B with concentration b on the 
left x G (— oo, — L) and on the right x G {L,oo) of the 
island. By symmetry our problem is reduced to the so- 
lution of the system (1) in the interval x € [0, oo) at the 
initial conditions 

a(x,0) = a o 0(L - x), b(x, 0) = b o 0(x - L) 

with the boundary conditions 

V(a,6)| x=0 =0, b(oo,t) = b . 

To simplify the problem essentially we will assume, as 
usually, Da = Db = D. Then by measuring the length, 
time, and concentration in units of L, L 2 /D, and bo, 
respectively, i.e. assuming L = D = bo = 1, and defining 
the difference concentration s(x, t) = a(x, t) — b(x, t) and 
the ratio of initial concentrations ao = r we come from 
(1) to the simple diffusion equation 



ds/dt 



with the initial conditions 



V 2 s 



(5) 



s (x e [0, 1)) = r, s (x e (l,oo)) 
and the boundary conditions 

Vs| x=o =0, s(oo,i) = -l. 
The solution to Eq.(5) has the form 



s(x,t) = 



r + 1 



erf 



1 + x 
2s/i 



erf 



l-x 
2Vt 



1, (6) 



whence, according to condition s(xf,t) = [2,19] there 
immediately follows the equation defining the law of mo- 
tion of the reaction front center, Xf(t) 



erf ( 1 + X J ) +crf ' 1 



2Vi 



2Vi 



r+ 1 



(7) 



Let us assume that d > d c and the reaction rate constant 
k is sufficiently large [30] so that at times t ~ toR oc 
k^ 1 <C 1 the annihilation goes to the scaling Galfi-Racz 
regime, i.e. in the vicinity of Xf there forms a narrow 
reaction zone w/An « 1. At t <C 1 and |1 — Xf \ <C 1 
from (7) in accordance with [2] we find 



x f 



1 



Vt- 



where erf(c//2) = (r — l)/(r + 1) and, hence, at r < 1 
the island contracts, whereas at r > 1 the island expands. 
By virtue of the fact that the number of particles in the 
island (per unit of the initial boundary) is finite, the stage 
of its expansion always goes to the stage of its subsequent 



contraction to end in an instant of time t c , when the 
reaction front center approaches the origin of coordinates, 
x f (t c )=0: 



evi(l/2VT c ) = l/(r + 1). 



(8) 



According to (8) at r ^> 1 the "lifetime" of the island 
t c 3> 1 so the majority of the particles die at times t ^> 1, 
when the diffusive length exceeds appreciably the initial 
island size, A D ^> 1. The evolution of the island in such 
a large-i regime is of principal interest to us here, and its 
analysis is the main goal of the present report. 
In the limit r,t> 1 from (6) we find 



s(x, t) 



VTTt 



(9) 



where \ = (1 — x 2 /2t)/12t H , and, hence, the law of 

the front motion is 



Xf = 2Vt(l + e)Wl 2 [( r -±l){l-e)], 



(10) 



where e = l/12t + ■ ■ -. By assuming f > 1 > ton,, at 
the stage of the island expansion and at a considerable 
time period of its subsequent contraction, the duration of 
which will be estimated below, the reaction zone width 
is w <C Xf < Ad, therefore in terms of the QSA we 
can take w/xf — ► 0, and, hence, R(x,t) = J5(x — Xf) 
so that a = s, b = when x < Xf, and a = 0, b = \s\ 
when x > Xf. Thus, by neglecting the width of the 
reaction zone, for the number of particles in the island we 
have N(t) = J Q Xf a(x,t)dx = f* f s(x,t)dx whence, with 
account taken of (9) and neglecting the contribution of 
X, we find 



N(t) = (r + l)ert(x f /2y/i) - x f . 



(11) 



Let us now calculate the maximal amplitude x 1 ^ 1 of 
the island expansion and the time tu of its achievement. 
According to (7) for r > 1 from condition if = we have 
exactly 



r¥ 
x f 



2t M Arcoth(a;f ) 



whence for r ^ 1 it follows 



(12) 



xf = V2^4l + l/12t M + ---)- 
Neglecting the terms 0(l/r 2 ) from (10), (12) we find 

t M = {r+ l) 2 /ne, xf = (r + \)^2jVe. 
Comparing tu with the lifetime of the island (8), 
t c - (r + 1) 2 /tt, 



we conclude that independently of r the ratio t c /tM — 
e = const. Going with allowance for this to the reduced 
coordinate £ = x/xf and time r = t/t c we, finally, come 
to scaling relations for the distribution of particles s((, r) 



2 



s (C,r)= e -< 2 / 2e 7V^-i, (13) 

for the center of the reaction front 

( f = v/er|ln<r|, (14) 
and for the number of particles in the island 

N/N = lr g(r), (15) 

where -f r = (r + — ► 1 as r — > oo and scaling function 

0(t) = crf(v/|lnr|/2) - v/2r| ln-r|/7r. 

From (13) and (15) we find that in the turning point 
t m = 1/e 

a(0,r M ) = v/e - 1 = 0.64872..., 

I^Nm/Nq = 0.19886..., 

and, hence, independently of the initial number of parti- 
cles, N = r, m 4/5 of the particles die at the stage of the 
island expansion and the remaining w 1/5 at the stage 
of its subsequent contraction. In Figs.l and 2 are shown 
the calculated according to Eqs.(6),(7) plots of the be- 
havior of the particle distribution |s((,t)|, of the center 
of the reaction front C/( t )j an d the number of particles 
r y~ 1 N((f)/N , which give the full insight into the scaling 
regime of the death of the island (13)-(15). 

Let us come now to the question of the scaling be- 
havior of the reaction rate R(x, t) in the vicinity of the 
reaction front center Xf. Assuming according to [8,9,19] 
that within the applicability of the inequalities tf/tj <C 1 
and w -C xj the reaction rate R(x,t) can be described 
in terms of the QSA by Eqs.(3),(4) to calculate Rf and 
w it only remains for us to calculate the diffusive current 
J(t) of particles arriving at the reaction zone. According 
to (9), (11) we have 

J = —N = —ds/dx \ X=XJ — Xf/2t 

whence it follows 

J/J M = Vllnrl/er, (16) 

where J M = 1/xf = y/ne/2/(r + 1). Substituting (16) 
into Eqs.(4) we come to the scaling of the reaction zone 

R f ~ [J^kf/* = Rf (^) 2/3 (17) 

w~ (Jk)-V 3 = w M (j^) , (18) 

where Rf ~ (fc/r 4 ) 1 / 3 and w M ~ (r/fc) 1/3 so that 
wRf <~ Rgiobai — Jg 00 Rdx — J. From (18) it is seen 



that at the stage of the island contraction the logarith- 
mic term becomes dominant, and as r — > 1, i.e., St = 
1 — t — > 0, the width of the reaction front w diverges as 
w cx (<5r) -1 / 6 . Comparing the characteristic times tj — 
-{dlnJ/dt)- 1 - r 2 5r and t f - w 2 - (r/fc) 2 / 3 ^)" 1 / 3 
weh&vetf /tj ~ (r 2 fc) _2 / 3 (<5r) -4 / 3 whence it follows that 
the characteristic time at which the QSA is violated is 

Stq oc 1/ry/k. 

Comparing then w and Xf we obtain w/xf ~ \/tf~Jtj ~ 
(r 2 fc) _1 / 3 ((5T) _2//3 whence it follows that the width of the 
reaction zone becomes comparable with the island size at 
the times 

5t w <~ 5tq cx 1/rVk. (19) 

From (19) it is seen that at sufficiently large values of 
i; » 1 the characteristic times 5t w ~ Stq <C l/r <C 1 
and, hence, practically all the particles die in the scaling 
regime (15). 

The analogous estimations in the case of the k- 
independent fluctuation Id regime with w ~ 1/VJ yield 
the expressions (17), (18) with Rf F - r~ 3 / 2 ,wfj ~ y/r, 
and exponents 3/4 instead of 2/3, and 1/4 instead of 1/6, 
respectively. In the end for the Id case we have 

8tZ ~ h$ cx l/r 2 / 3 . (20) 

According to (20), despite the much more "blurred" 
structure of the Id island, in this case, too, in the limit 
of sufficiently large r » 1 the majority of the particles 
die in the scaling regime (15). 

In summary, the problem of death of an A particle is- 
land in the B particle sea at equal diffusivities of A's and 
B's particles has been first considered. It has been found 
that at sufficiently large initial number of particles in the 
island, r » 1, and sufficiently large reaction rate con- 
stant, k 3> 1, the death of the majority of the particles at 
any d occurs in the universal scaling regime, and the most 
essential features of this regime have been revealed. The 
obtained results can have many applications, especially in 
surface science, and, in particular, they can provide fresh 
insight into the dynamics of the Ovchinnikov-Zeldovich 
hierarchic A — B structures [1]. By analogy with [19] it 
may be expected that the analysis presented can be eas- 
ily extended over the general case Da 7^ f b ^ (note 
that in the static case Db — 0, which belongs to the sep- 
arate universality class, with initially N particles A at 
the origin of coordinates, in [16] the formally coincident 
with (15) scaling, N = N^lt/N^), was found). An in- 
vestigation into the general case of nonzero diffusivities 
and a most interesting generalization for the anisotropic 
diffusion of reactants is expected to be presented in a 
future report. 

I would like to thank Alexey Nekrasov for help in com- 
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the RFBR through Grant No. 02-03-33122. 
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FIG. 2. (a) Inset: Time dependences Xf(t) calculated from 
Eq.(7) at r = 5 (filled squares), r = 10 (open squares), r = 20 
(filled circles), and r = 40 (open circles). Main panel: Col- 
lapse with the growing r of the shown in the inset dependences 
to the scaling law (14) (solid line) in the rescaled coordinates 
(/(t). For completeness is shown the motion of the both 
island fronts: ±Xf(t), ±C/(r); (b) Collapse with the grow- 
ing r of the calculated from Eqs.(6) and (7) dependences of 
the reduced particle number in the island 7 r r 1 jV/7V vs to 
the scaling function <J(£/)(solid line): r = 5 (filled squares), 
r = 10 (open squares), r = 20 (filled circles), and r = 40 
(open circles) . The number of the island particles N has been 
calculated by integrating Eq.(6) from to i/ on assuming 
that w/xf — » 0. 



FIG. 1. Evolution of the distribution of particles |s(£, t)| 
calculated according to Eqs.(6) and (7) for r = 10. The left 
part of the distribution (£ < 0) is obtained through the mir- 
ror reflection of the right one (£ > 0). The region s > 
(— C/ < C < C/)i belonging to the island, is colored grey. For 
demonstration only a part of the picture is shown for r > 0.03. 
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